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Using ah initio methods based on the density functional theory, we investigate the segregation 
and formation energies for various dopants (Si, Be, Zn, Sn), commonly used to obtain p- or n- 
type conductivity in GaAs and InAs nanowires. The distribution of Au and O atoms, which may be 
unintentionally incorporated during the wire growth, is also studied. The calculations performed for 
nanowires of zinc blende and wurtzite structure show that the distribution of most of the impurities 
depends on the crystal structure of the wires. For example, it is shown that the same growth 
conditions can lead to lower energy for Si substituting Ga (donor) in the wire of wurtzite structure 
and substituting As (acceptor) in the wire of zinc blende structure. In contrast, we obtain that 
gold and oxygen atoms always tend to stay at the lateral surfaces of GaAs and InAs nanowires, 
in agreement with experimental findings, while for beryllium the lowest energies are found when 
the impurities are located in sites in the center of the wurtzite wire or along the [1, 0, 1] axis from 
surface to the center of the zinc blende wire, what can explain the recently observed diffusion of this 
impurity into the volume of GaAs wires. 

PACS numbers: 63.20dk, 61.72uj, 62.23Hj, 61.72sh 



I. INTRODUCTION 

Successful realization of most of electronic nanodevices 
requires a high carrier concentration so as to reach suf- 
ficient conductivity. The excellent properties of III-V 
semiconductors, in particular the high electron mobility 
which can be obtained in GaAs and InAs, turns nano- 
wires (NWs) from these materials into natural candi- 
dates for building blocks of such high speed applications. 
The properties of semiconductor NWs diflFer from those 
of their parent bulk materials due to quantum confine- 
ment and surface related effects - pure semiconductor 
NWs usually have low concentration of carriers. Thus, 
the application of semiconductor NWs in novel electronic 
devices requires controllable p-type and n-type doping. 
This task still remains a big challenge to the grow- 
ers, despite substantial eflForts, both experimental and 
theoretical.— ^ Most of the publications describing dop- 
ing of III-V NWs focus on InAs NWs.^ In these structures 
the surface Fermi level is pinned in the conduction band. 
This makes n-type conductivity in InAs NWs easy to ob- 
tain at the expense of difficulties in p-type doping.— In 
contrast, in GaAs NWs the Fermi level at the surface is 
pinned approximately in the center of the band gap and 
both, controlled p- and n-type doping, can be produces. 
Doping of GaAs NWs grown by molecular beam epitaxy 
(MBE) has been demonstrated in different means. Cz- 
aban et al. used Be and Te as p- and n-type dopant 
precursors respectively^ while Fontcuberta i Morral et 
al. pointed out that Si may act as both, acceptor and 
donor, by just changing the operating temperature dur- 
ing growth as shown before for 2D epitaxial growth of 
GaAs.^^'^^ The incorporation of Si and Be into GaAs 
NWs was investigated in a further study Recently, the 
use of tetraethyl tin as a precursor material for successful 
n-type doping of GaAs NWs, grown by metal organic va- 
por phase epitaxy, was reportedJ^ Moreover, it occurred 
that not only specific growth conditions but also the crys- 



tal structure of the wires may lead to different incorpo- 
ration behavior than known for planar layer sJ^— This 
was shown in particular for the group IV element, silicon, 
which is known to be amphoteric in III-V compounds. 
Namely, it was found that although under typical MBE 
growth conditions Si acted as n-type dopant for GaAs, 
p-type conductivity was obtained in Si-doped GaAs NWs 
which adapted wurtzite (wz) structure in similar growth 
conditions i 

As mentioned above, doping of NWs was studied the- 
oretically by several groups. They focused, however, 
mainly on silicon NWs (see Ref.[l9| and the references 
therein). There are very few publications which describe 
theoretically, incorporation of non-magnetic dopants in 
GaAs and InAs NWs. To our best knowledge, doping of 
GaAs NWs was considered only in Ref. [H, where only Si 
doping of wz structure GaAs NWs was studied. In con- 
trast, Claudia dos Santos and co-authors considered Cd 
and Zn doping of InAs NWs having the zinc-blende (zb) 
structure.^ In Ref. the studied NWs were extremely 
thin (with diameters of 1 nm) and their lateral surfaces 
were passivated with hydrogen atoms. Recently, the role 
of the surface danghng bonds and the effect of molecular 
passivation on doping of InAs NWs was also studied by 
Haibo Shu et alf^^ 

In this paper we report a systematic theoretical study 
of several dopants, which are usually used to obtain p- 
or n-type conductivity in GaAs and InAs NWs, i.e., Si, 
Be, Sn and Zn. In contrast to the previous papers, we 
assume that the surface atoms are not saturated by fo- 
reign atoms. This assumption is justified by the fact, that 
we consider incorporation of dopants during the growth 
process in high vacuum conditions, e.g., MBE. To check 
what is the preferential location of a given dopant in the 
wire, i.e., to answer the question whether it is possible to 
explain the problems of doping of III-V NWs by segrega- 
tion phenomena, we study the energy of the wire with the 
impurity situated in diflFerent sites. Such analysis proved 
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already to be very useful, e.g., for predicting conditions 
to obtain ferromametism in GaAs NWs with magnetic 
Mn ions^^ In Ref. [20| it was shown that the distribution 
of Mn ions and thus the electric and magnetic proper- 
ties of (Ga,Mn)As NWs depend strongly on the crystal 
structure. Since the III-V semiconductor NWs can grow 
in both, zb and wz structures, we consider here both 
types of wires to check whether the crystal structure has 
an impact on the distribution of the studied impurities 
in GaAs and In As NWs. Finally, we consider the distri- 
bution of Au and O atoms in the studied NWs. Oxygen 
can be used for n-type doping of GaAs— but also un- 
intentional incorporation of O atoms in GaAs wires has 
been reported in the literatures^ On the other hand, gold 
can defuse into the wire from the catalyst particle during 
Au-assisted vapor - liquid - solid (VLS) growth. This is 
known to be a serious problem in Si NWs.^^ Although 
it seems that the unintentional incorporation of Au in 
GaAs and InAs wires is much lower^ the measured Au 
doping values can make ballistic transport through the 
nanowires difficult to obtains^ This makes the study of 
the distribution of Au atoms in the wires very important. 
It should be mentioned here that to avoid Au incorpo- 
ration self-catalyzed growth of GaAs and InAs NWs has 
been pursued by several groups (see, e.g.. [27I - [3q1 ). 

II. METHODOLOGY 

In our calculations the doping of the wires is performed 
by substituting one cation or anion in the elementary cell 
of the NW with an impurity atom. One foreign atom in 
the elementary cell of the studied NWs corresponds to a 
concentration of impurities in the range from 0.5% to 
^1.6%, depending on the crystal structure and diame- 
ter of the wire. We study the distribution of impurities 
by comparing the total energies of the NWs with an im- 
purity at various non-equivalent substitutional sites. In 
Fig. [T] we present the dopant positions considered here. 
To a given cation site and the anion above this cation we 
ascribe the same number. It should be emphasized that 
in the cross section corners, i.e., in the positions 1 and 6, 
the geometry of bonds is different for cations and anions. 
Namely, in site 1 the cation has two dangling bonds while 
the anion has one dangling bond, and vice versa in site 
6. 

The energies of the NWs are calculated by ab initio 
methods, based on the density functional theory (DFT). 
The energy of the wire with the dopant occupying the site 
exactly at the center of the wire cross section is taken as 
a reference point, i.e., as zero on the energy scale. The 
segregation energy is defined by the difference between 
the energy of the NW with an impurity at a given site 
and the reference energy: 

= E doped NW — E center doped NW - (1) 

Thus, the negative value of the segregation energy for 
a dopant atom at a given site means that this impurity 




FIG. 1. (Color online) Cross section of an a) wz NW with 
diameter d=2.8 nm and b) zb NW (d=2.4 nm). The green 
balls show the 10 nonequivalent sites, which we consider as 
substituted by a dopant. The blue balls denote Ga(In) atoms, 
gray - As atoms. The darker colors denote atoms in the deeper 
double layer. 



prefers to substitute the host atom at this particular site 
rather than be at the center of the wire. In contrast, 
positive segregation energy values denote the situation 
when substitution of an atom in the center of the NW 
will be energetically more favorable. Let us mention that 
in a study of doped Ge wires, performed with use of 
a similar method, the values of formation energy, i.e., 
energy needed to insert the impurity atom (taken from a 
reservoir) into a given NWs site, from which the original 
atom has been removed (to a reservoir), were shown to 
depend strongly on the size of the simulated supercell and 
that segregation energy, being a relative quantity, is free 
of these systematic inaccuracies. However, by the study 
of segregation energies we can only answer the question 
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where an impurity prefers to stay within the wire. To 
answer the question for which impurity the energy cost 
of substituting a host atom is lower, one has to compare 
their formation energies, which depend on atomic and 
electronic chemical potentials. The formation energy Q 
of a neutral impurity X can be calculated according to 
the formalism proposed by Northrup and Zhang^ 

n = Etot [X] - Etot [NW] + ^Ga/As - MX , (2) 

where Etot [X] is the total energy derived from the calcu- 
lation of the nanowire with an impurity X, Etot[NW] is 
the total energy of undoped nanowire. The term fica/As 
denotes the chemical potential of the corresponding sub- 
stituted Ga or As atom, while /ix denotes the chem- 
ical potential of the impurity atom. Our analysis is 
performed with an assumption of thermal equilibrium 
conditions of the impurity region with the surrounding 
nanowire, i.e., 

MGaAs = MCa + MAs. (3) 

Using fii < with i G Ga, As, X and the definition 

of the heat of formation of GaAs 

A U ,,bulk I .hulk .hulk / a\ 

we can limit the range of possible values of A /i to — AH < 
A/i < -\-AH. The quantity A/i is the difference between 
the chemical potentials of Ga and As and is given by 

Am = i^lGa - has) - (A^^-f - hTs')- (5) 

Thereby, the situation A/i = -\-AH corresponds to Ga- 
rich growth conditions, while A/i = —AH corresponds 
to As-rich conditions. 

As mentioned above, we incorporate the dopants into 
GaAs and InAs wires of both wz and zb crystal struc- 
tures. The only considered NWs are those grown along 
[0001] direction (wz) and along [111] axis (zb), because 
a vast majority of III-V NWs grow only along these di- 
rections. It has been also shown theoretically that these 
growth directions and the types of facets presented in 
Fig. [1] constitute the most stable configurations of GaAs 
and InAs NWs*^ — It is well known that with ab initio 
methods only very thin NWs, much thinner than the real 
ones, can be simulated. This can lead to overestimation 
of the side surfaces' effects in the results, deceptive when 
such calculations are used to determine, e.g., the crys- 
tal structure of thicker wires. This problem is, however, 
not so important in our study, in which we try to answer 
only the question: is there any energy difference between 
the wires with an impurity situated at the side facets and 
with the impurity inside the core of the wire. The diame- 
ter test, in which we considered zb wires of several diffe- 
rent diameters, 0.9 nm, 1.4 nm, 1.8 nm and 2.4 nm, has 
shown that for the two highest diameters the obtained 
segregation energies are practically the same. Similarly, 
the results do not depend considerably on the diameter 
in wz NWs with diameters from 2.0 nm to 2.8 nm. Con- 
gruous studies reported in Ref. [H also indicate that for 



diameters larger than approximately 3 nm the P impurity 
levels inside the Si NWs recover bulk characteristics. In 
the following, all the presented results are for the thick- 
est wz NWs with diameters of 2.8 nm and the thickest 
(2.4 nm) zb wires. As already mentioned, we assume 
that the surface atoms are not saturated by any foreign, 
for example, hydrogen atoms, because we study the pro- 
cess of NWs growth and incorporation of the impurities. 
Still, for the calculations of the density of states of al- 
ready doped wires and discussion about the positions of 
the Fermi level, the side facets of the wires have been 
passivated by partial atomic charges of hydrogen atoms. 

The calculations are performed using the Vienna 
ab initio simulation package (VASP)*^i^ Each of the 
initial structures is cut out from an appropriate bulk ma- 
terial, in either zb or wz structure. The wz elementary 
cell is doubled along the c direction, to guarantee that 
there is enough separation between the impurity ion and 
its periodic images. Thus, the smallest supercell used in 
our study contains 122 atoms (zb) whereas the largest 
(wz) 384 atoms. For different dopants a plane wave basis 
set is adopted with different cutoff energies: for Si the 
cutoff energy was equal to 307 eV, for Be - 375 eV, Sn 
- 261 eV, Zn - 346 eV, Au - 287 eV and O - 500 eV. In 
each case all the energies and chemical potentials are cal- 
culated with these cutoff values. These cutoff energies are 
sufficient to ensure that the results do not depend on the 
size of the plane wave basis. The k-points are generated 
with a (1 X 1 X n) mesh, where n > 50 A/c and c is the 
unit cell dimension in the growth direction. The inter- 
action between the valence electrons and the ionic cores 
is included using the projected augmented wave method 
(PAW)*^ The exchange correlation energy is calculated 
within the generalized gradient approximation (GGA). 
The atomic positions leading to minimum energy are de- 
termined with relaxation of all atomic positions as well 
as of the unit cell dimension in the growth direction, and 
with full reconstruction of the NW's surfaces. In each 
simulation process the atomic coordinates are relaxed 
with a conjugate gradient technique. The criterion that 
the maximum force is smaller than 0.01 eV/A is used to 
determine the equilibrium configurations. Neither pres- 
sure nor temperature are included in our calculations. 



III. RESULTS 

The distribution of gold in the cross-section of InAs 
wires was studied experimentally using a technique called 
pulsed-laser atom probe tomography. In Ref.|25| it was 
shown that Au atoms accumulate in a shell at the lateral 
surfaces of the wire. We start our study by considering 
Au impurities in InAs NWs to check if our calculations 
can describe this experimental finding. We analyze InAs 
NWs with an Au atom in various cation sites and focus 
on the search of the impurity positions for which the 
segregation energy is the most negative, i.e., the sites in 
which the energy cost of substituting the cation by an 
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Au atom is lower than in the other nonequivalent sites. 
In Fig. [2] the results for both wz and zb InAs as well 
as GaAs NWs are presented. In this and all following 
Figures the diameters of the wz NWs equal to 2.8 nm 
and of zb NWs to 2.4 nm. In all Figures the calculated 
values of segregation energy are connected by a line to 
guide the eye. 
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FIG. 2. (Color online) The segregation energy in a) wz (dia- 
meter 2.8 nm); b) zb (diameter 2.4 nm) InAs (green) and 
GaAs (blue) NWs as a function of different Au atom posi- 
tions. 

As one can see in Fig. [21 in both wz and zb NWs the 
Au atoms tend to accumulate near the lateral surfaces. 
In wz structure the energy is the lowest for the dopant 
occupying the threefold coordinated sites 1 and 4, the 
corner sites of the sidewall facets. In NWs of zb struc- 
ture the most energetically favorable for Au atom is to 
substitute the cation in position 1, i.e., the cation with 
extra dangling bond in the corner of the cross-section of 
the wire. Analogous results are obtained for oxygen, an- 
other element, which can be unintentionally incorporated 
during the growth, as shown for GaAs wires in Fig. [3l 
The difference between the energy of the wire with an 
Au or O atom at the lateral surface and at the center 
of the wire is significant, from ca 0.6 eV for gold in wz 
InAs NWs to nearly 3 eV for oxygen in GaAs zb NWs. 
Thus, both these impurities should indeed be trapped at 
the lateral surfaces of the wires. In zb NWs we note that 
the segregation energy decreases with the distance from 
the center of the wire monotonically along the different 



crystallographic axes, as shown for example in the inset 
to Fig. [3] for two nonequivalent [1,0,1] directions. 
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FIG. 3. (Color online) The segregation energy in a) wz; b) 
zb GaAs NW as a function of different O atom positions. In 
the inset to b) the segregation energies for the O atom in 
subsequent sites along the two nonequivalent [1, 0, 1] lines are 
shown 

Next, we consider Si doping in GaAs NWs. As already 
mentioned, since Si is a group IV element in III-V com- 
pounds it can serve either as a donor, by substituting 
the cation, or as an acceptor, by occupying the anion 
site. We therefore analyze the NWs with the Si atom 
in various cation and anion sites. In Fig. |4] the results 
for both wz and zb NWs are presented. We recall that 
Si-doped GaAs NWs of wz structure were already stu- 
died by Ghaderi et al. in Ref. 0. These authors con- 
sidered only four dopant positions in a wz NW, i.e., one 
center-like, two surface threefold (distinction is made for 
corner or middle of facets positions), and one subsurface 
fourfold coordinated position. Our results presented in 
Fig. 13^, which are obtained for many more Si positions 
in a thicker wz GaAs NW, fully confirm the conclusion 
of Ref. 8. Namely, the four-coordinated subsurface po- 
sitions are the most energetically favorable locations for 
Si acting as a donor as well as an acceptor. The seg- 
regation in these sites is, however, not as strong as for 
Au or O atoms (ca -0.2 eV). The calculated segregation 
energies for Si-doped GaAs NWs of zb structure, which, 
we recall, were not considered by the authors of Ref. [l, 
are presented in Fig. HJ). As one can notice by compar- 
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FIG. 4. (Color online) The segregation energy in a) wz; b) zb 
GaAs NW as a function of different Si atoms positions. 



ing Fig. EH and b, the spread of the segregation energy 
values is considerably larger for zb wires (±0.7 eV). It 
means that one can expect much more homogeneous dis- 
tribution of Si atoms in wz wires than in those of zb 
structure. For both Si donors and acceptors the most 
preferable position in zb GaAs wires is site 7, i.e., the 
fourfold coordinated site next to site 1. However, for Si 
substituting Ga (donor) low energy is also obtained for 
site 11, in the vicinity of the corner anion with an extra 
dangling bond. We can thus conclude that our calcu- 
lations suggest that while in wz GaAs NWs the energy 
barrier at the lateral surface is not very high and Si can 
be incorporated more or less homogeneously, in zb GaAs 
NWs Si doping should result in accumulation of Si at the 
shell. 

The impact of the crystal structure on the incorpo- 
ration of Si into the III-V NWs is presented in Fig. [5l 
In this figure the formation energies for Si substituting 
Ga and As in the center (site 16) of the wz GaAs NW 
(Fig. El a) and center (site 21) of zb GaAs NW (Fig. [5] 
b) are compared. We observe that there exists a region 
of A/i in which the same growth conditions lead to lower 
energy for Si substituting Ga in the wire of wz structure 
and substituting As in the wire of zb structure. 

Next, we have calculated the density of states (DOS) 
for Si-doped GaAs nanowires. In wz structure wires Si 
acts as expected, i.e., as a donor when substituting Ga 
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FIG. 5. (Color online) The formation energy Q of GaAs NW 
with Si impurity in the center of a) wz NW, b) zb NW as a 
function of the value A/i. 



and as acceptor when substituting As. Impurity bands 
are located close to valence band maximum (VBM) and 
conduction band minimum (CBM) for Si acceptors and 
donors, respectively. The passivation of the lateral sur- 
face have almost no influence on the position of these 
states, similarly to what was shown for CdTe NWs of wz 
structure in the paper of T. Sadowski4S In contrast, in 
not passivated zb NWs we observed surface states in the 
band gap and the impurity states (for both Si at As and 
Ga sites) located at these surface states. After satura- 
tion the surface states disappear and the impurity states 
move close to VBM and CBM, for Si atom substituting an 
anion or cation atom, respectively, inside the wire core. 
In Fig. [6]a) we show the density of states calculated for 
H-passivated zb GaAs nanowire with Si substituting the 
Ga atom in the wire's center position 21. The situa- 
tion is quite different when Si substitutes the Ga atom in 
position 1, in the corner of the wire cross-section. For un- 
passivated nanowires (part b of Fig. [6]) Si states and the 
Fermi level are again pinned at the surface states in the 
band gap, as for the other impurity positions in zb NWs. 
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FIG. 6. (Color online) The density of states calculated for 
zb GaAs NW with Si substituting Ga atom in the site 21 a) 
and site 1, c) and b), (site numbers as denoted in Fig|T]). The 
DOS presented in part b) is for the NW with unpassivated 
lateral surfaces; the results in parts a) and c) are for NWs 
with lateral surfaces passivated by H atoms. Zero on energy 
scale denotes always the Fermi level. 



However, after passivation the Fermi level moves to the 
VBM (compare part c of Fig. [6|), what indicates that 
in this position Si will not behave as an n-type dopant. 
As the segregation energies suggest that most of the Si 
atoms are trapped around these corner positions, this be- 
havior can explain the observed dominance of p-type in 
Si-doped GaAs NWs.J^ 

In the next step we compare the segregation energy 



obtained for Si-doped GaAs NWs with the segregation 
energy in wires doped with beryllium, one of the most 
widely used acceptors in GaAs. Results obtained for Be- 
doped GaAs NW of wz and zb structure are presented in 
Fig. [71 As shown in the Figure, in wz NWs the distribu- 
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FIG. 7. (Color online) The segregation energy in a) wz; b) zb 
GaAs NW as a function of different Be atom positions. 

tion of Be atoms should be again quite homogenous, but 
in this case the lowest energy corresponds to the situa- 
tion where the impurity substitutes an atom inside the 
wire. Again, in zb GaAs NWs the lowest value of the 
segregation energy (-0.4 eV) was obtained when Be sub- 
stitutes the Ga atom with an extra dangling bond in site 
1. It should be noted, however, that this is much weaker 
segregation than obtained in zb GaAs NWs for Si, and of 
course Au and O. Therefore, these results seem to be in 
agreement with the experimental observation that even 
though Be atoms accumulate in the shell of zb GaAs 
NWs, it should be possible to dope the NWs almost ho- 
mogeneously during axial growth. It is also interesting 
to note that for some reason the lowest energies of the zb 
GaAs wires with Be atoms are found when the impuri- 
ties are located in sites along the [1,0, T] axis (from the 
center of the wire to the cation with extra dangling bond 
at the corner of the lateral surface). This is in contrast 
to the situation presented in Fig. [3] for oxygen and gold 
impurities. The low energy cost of incorporation of Be 
along the [1,0,1] axis can explain the recently observed 
diffusion of this impurity into the volume of the zb GaAs 
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NWs, interpreted in Ref. '4l' as an alternative incorpora- 
tion path. 

Finally, the comparison of formation energies for Be 
incorporation into the GaAs NWs, presented in Fig. [51 
shows that it should be much easier to dope with beryl- 
lium the NWs with zb structure than those with wz crys- 
tal structure. It shows also that the formation energies 
for Be are much lower than for Si. 
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FIG. 8. (Color online) The formation energy Q of zb and wz 
GaAs NW with Be impurity as a function of the value Afi. 

As mentioned before, the as grown In As NWs are usu- 
ally n-type. Doping with several acceptors, like Be or Zn, 
have been tried to obtain p-type conductivity in these 
NWsi^i^ In Fig. [9] we present the comparison of segrega- 
tion energies obtained for InAs NWs doped with Si, Be 
and Zn atoms as a function of their different locations. 
The study of segregation energy for different acceptors in 
InAs NWs shows that in these wires, similarly to GaAs 
NWs, wz structure should lead to a more homogeneous 
distribution of impurities. Inside the wire, in the core, 
all impurity positions are energetically equivalent in both 
structures. In the shell. Be and Zn prefer to be in the 
middle of the sidewall of the wire, substituting three- 
fold coordinated cations, while Si, substitutes fourfold 
coordinated anion at the NW's lateral surface (site 6), 
as presented in Fig. [TOl One can understand this dif- 
ference by recalling that Zn and Be occupy the cation 
positions while Si, should substitute an anion to serve 
as an acceptor. The relaxation of the atom positions at 
the wire lateral surfaces leads to different reconstruction 
of bonds for cations and anions, as shown in Ref. [sH for 
pure III-V NWs. The most important result of this study 
is, however, the fact that in InAs NWs of wz structure 
the segregation energy for Be atoms is not negative for 
any of the considered positions. Also for Si in wz InAs 
NWs the segregation energy for the impurity in site 6, 
although negative, is close to zero (-0.06 eV). Thus, we 
can conclude that the best precursor for p-type doping of 
InAs NWs should be beryllium but using Si can also be 
effective, provided the wz crystal structure of the wires 
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FIG. 9. (Color online) The segregation energy for Si, Be and 
Zn acceptors in various positions in wz (a) and zb (b) InAs 
NW. 




FIG. 10. (Color online) The most energetically favorable po- 
sitions of Zn, Be and Si acceptors at the lateral surfaces of wz 
InAs NW. 



is assured. 

In Ref. Hslit was shown, by capacitance measurements, 
that the concentration of carriers in InAs NWs can be 
increased by n-type doping with Sn. In Fig. [11] the seg- 
regation energy obtained for InAs NWs doped with Sn is 
presented. The results for Sn in wz InAs wires are very 
similar to that for Si-doped InAs NWs, presented above. 
Tin ions distribute fairly homogeneously across the wz 
wire. Although for Sn the most favorable is to substitute 
the fourfold coordinated cations at the wz NW's late- 
ral surface, the segregation energy is again very small. 
This result seems to agree with the observed increase of 
the carrier concentration and the surface charge density 
in Sn-doped wz InAs NWs.^^ It should be noted that 
in Ref. 3 it was shown that the Sn precursor increases 
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FIG. 11. (Color online) The segregation energy of a) wz; b) FIG. 12. (Color online) The segregation energy of a) wz; b) 



zb InAs NW as a function of different Sn positions. 



zb GaAs NW as a function of different Sn atoms positions. 



the stacking fault density in wz NWs, ultimately at high 
flows leading to a zb crystal structure with strong over- 
growth and very low resistivity. In agreement with this 
observation we obtain that in zb NWs, analogously to 
other dopants in this material, Sn atoms do not incor- 
porate into the wire but stay at the lateral surface, with 
particular reference to atoms with extra dangling bond 
(Fig.IIIb). 

In contrast, our calculations performed for GaAs NWs 
show that in this material Sn atoms tend to segregate 
to the lateral surfaces in both wz and zb structures (see 
Fig. [12]). Thus, although our calculations performed for 
Sn impurities in III-V wires agree with the observation 
that doping with tin can be used to increase the con- 
centration of electrons in InAs wires^ they suggest that 
using Sn for n-type doping of GaAs NWs can be much 
less efficient. 



IV. CONCLUSIONS 

Nanowire growth processes usually rely on the vapor - 
Hquid - solid (VLS) mechanism normally assisted by gold 
droplets or self assisted on a SiO/Si02 surface. Doping 
these structures during such complex growth process is 
difficult to understand and control. It is, however, al- 
ready widely accepted that there are two distinct path- 



ways for dopant incorporation into semiconductor NWs: 
the VLS mechanism, in which the dopant atoms first dis- 
solve in the catalyst and are then transported across the 
hquid - solid interface to the core of the NW, and the 
vapor - solid (VS) mechanism, when dopant atoms are 
directly deposited on the lateral NW surfaces. In Ref. 
Il44 it has been suggested that these two mechanisms pro- 
ceed at different rates - namely, the latter predominates 
and leads to a much higher concentration of dopants in 
the shell of the NW than in the core. This conclusion 
was based on the experimentally observed accumulation 
of gold atoms close to the lateral surfaces of InAs NWs— 
as well as phosphorus in Ge NWs.^^ Radovanovic and co- 
workers, when studying GaN wires doped with Co and 
Cr, have also concluded that adsorption on the lateral 
surface is the main doping mechanism of the semicon- 
ductor NWs 4^ Finally, the conclusion that Be atoms are 
mostly incorporated into the zb GaAs NWs from the side 
facets and that the incorporation through the Ga droplet 
is negligible are made in Ref. |4l|. It should be also men- 
tioned that the very recently reported beautiful model^ 
in which B-doped Si NWs were considered, has shown 
that 2D diffusivity along the solid-liquid NW interface 
plays also a considerable role in impurity incorporation. 
In this paper we show, by means of first-principles cal- 
culations of the total energy of doped GaAs and InAs 
NWs, that indeed most of the commonly used impurities 
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for p- and n-type doping should accumulate near the sur- 
faces of the NW. It should be, however, emphasized that 
this result does not depend on the growth mechanism, 
because in our calculations only the differences in for- 
mation energies for various dopant positions have been 
taken into account. Thus, it shows that even for incorpo- 
ration of the dopants into the core, via VLS mechanism, 
the segregation of impurities to the surfaces should take 
place. The only exceptions are wz GaAs as well as wz 
InAs NWs doped with Be - in these the segregation en- 
ergy is positive for all studied cation sites substituted by 
Be atoms. 

We note that the degree of segregation depends con- 
siderably on the NWs host material, crystal structure of 
the wire and of course the particular dopant. Our calcu- 
lations show that in NWs of zb structure the segregation 
for the most common dopants is much higher than in 
NWs of wz structure. The studied impurities in zb NWs 
should remain either at the lateral surfaces, where they 
prefer to substitute the atom with extra dangling bond, 
or should be trapped at the subsurface, in the vicinity of 
the atom with additional dangling bond. At these posi- 
tions the segregation energies in zb GaAs NWs are lower 
than in the core by ca 0.4 eV for beryllium, but by even 
1.5-3 eV for gold and oxygen. It is also shown that at 
these positions Si atoms do not act as donors even when 
substituting Ga cations - this result is in agreement with 
the observed predominance of p-type in Si-doped GaAs 
NWs. For gold and oxygen in wz wires the results are 
similar to those obtained for zb structure but with lower, 
though still high, segregation. With our calculations we 
have thus confirmed the experimentally observed segre- 
gation of Au atoms to lateral surfaces^^ According to 
our results the same should be also observed for oxy- 
gen. In contrast, all other studied impurities distribute 
much more homogeneously across the wz wire and Be 
even prefers to substitute the Ga(In) ion at the center of 
the wire, as mentioned above. We can thus conclude that 
our calculations suggest that growth conditions leading 
to wz structure can help to avoid accumulation of im- 
purity atoms at the surface during the growth of doped 
III-V wires. It seems that the best choice for effective 
p-type doping of both GaAs and InAs wz NWs is to use 
berylhum as a dopant, due to it's low segregation. How- 
ever, the cost of substituting the host cation by Be atom, 
i.e., the formation energy, is lower in the zb than in the wz 



NWs. Finally, we observe that the different from other 
impurities behavior of Au and O can be attributed to the 
large electronegativity of these elements (3.5 for O and 

2.4 for Au) as compared to electronegativities of Ga, In 
and all other studied dopants, which are in the range of 

1.5 - 1.8. 

As shown above, our results suggest that the energy of 
a doped NW can be very different when the impurity po- 
sition inside the wire changes, especially in zb structure. 
Only in the case of oxygen and beryllium in GaAs NWs 
some regularities in these differences can be observed. 
In other cases the fluctuations of the segregation energy 
values seem to be random. These results were fully con- 
firmed by several tests, i.e., by repeating the calculations 
for the impurities in equivalent sites, by changing the 
NW diameter and finally by moving slightly the initial 
atomic positions. Recently, thanks to laser assisted atom 
probe tomography the studies of the distribution of impu- 
rities in semiconductor wires have become possible. Our 
results may suggest that the observed inhomogeneities 
(compare, e.g., Ref. [l9l ) may have a more fundamen- 
tal origin than just disorder. Neither the noticed trends 
nor the other segregation fluctuations, which we are not 
able to explain in our approach, were reported before. 
It should be noted, however, that in all previous studies 
the segregation and/or formation energies in doped NWs 
have been calculated for at most 3-4 different impurity 
positions. In our study the results for wz NWs with the 
impurities in 16 different sites and for zb NWs in 21 dif- 
ferent sites have been compared - to our best knowledge, 
most of the nonequivalent sites considered here have not 
been studied before for any impurity. 
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